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Abstract 

For materials science, diamond crystals are almost unrivaled for hard¬ 
ness and a range of other properties. Yet, when simply abstracting the 
carbon bonding structure as a geometric bar-and-joint periodic frame¬ 
work, it is far from rigid. We study the geometric deformations of this 
type of framework in arbitrary dimension d, with particular regard to the 
volume variation of a unit cell. 
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Introduction 

In this paper we survey a number of topics in the deformation theory of pe¬ 
riodic bar-and-joint frameworks by investigating diamond frameworks, that is, 
generalizations to arbitrary dimension d of the ideal atom-and-bond structure 
of diamond crystals. 

Our setting is purely geometric and does not involve any physical assumptions or 
properties. A general deformation theory of periodic frameworks was introdnced 
in [3] and further developed in HIE]- Although crystallography, solid state 
physics and materials science have studied crystal structures for a long time, 
abstract mathematical formulations are of relatively recent date m- Of course, 
there are historical roots in problems related to sphere packings [UllillH], 
lattice theory [3 [23] and crystallographic groups miEiiii]. 

Diamantine frameworks have elementary geometrical descriptions in any dimen¬ 
sion d. The standard planar case d = 2 is the familiar hexagonal honeycomb 
illustrated in Figure Another configuration, the so-called reentrant honey¬ 
comb illustrated in Figure]^ was recognized in materials science as a structure 
with negative Poisson’s ratio m [22] and has become an emblem of auxetic 
behaviour [niiiT]. It will be seen below that auxetic deformation paths can 
be defined in strictly geometric terms [7] and resemble causal-lines in special 
relativity. 

Our treatment is mostly self-contained and may serve as an “introduction by 
example” to several topics of general interest: topology of the deformation space. 
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variation of volume per unit cell and critical configurations, possibilities for 
auxetic trajectories. 



Figure 1: The standard diamond framework in dimension two. The generators 
of the periodicity lattice are emphasized by arrows. The quotient graph has two 
vertices connected by three edges. 



Figure 2: A reentrant honeycomb is a deformed diamond framework with 
auxetic capabilities. A horizontal streching detrmines a vertical expansion. 


1 The standard diamond framework in 

The standard or canonical diamond framework in R'^ can be described by start¬ 
ing with a regular simplex PqPi...P^ centered at the origin O. Then, we take 
the midpoint Mq of the segment OPq and denote by QoQi---Qd the simplex 
obtained from the original one by central symmetry with center Mq. Figure 
shows this setting in dimension three. By using as periodicity lattice A the rank 
d lattice generated by the vectors Xi = Pi — Pq, i = 1, ...,d, we can produce a 
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periodic set of vertices which includes the vertices Pi and Qj, by considering 
all translates of Pq and Qo = O. For edges we take all segments OPi — QoPi, 
i = 0, ...,d and their translates under A. Up to isometry and rescaling, the 
resulting framework is unique. 



Figure 3: An illustration for d = 3. The standard diamond framework uses 
the vertices of two regular tetrahedra. Each tetrahedron center is a vertex of 
the other tetrahedron and the pair is cetrallly symmetric with respect to the 
midpoint of the two centers. Centers are seen as joints connected by bars to the 
vertices of the other tetrahedron. By periodicity under the lattice of translations 
generated by all edge-vectors of the tetrahedra, one obtains the full bar-and-joint 
diamond framework. 

The abstract infinite graph together with the marked automorphism group 
represented by A gives a d-periodic graph as defined in Hi]. There are two 
orbits of vertices and d+1 orbits of edges under this periodicity group. In other 
words, the quotient multigraph consists of two vertices connected by d-|-1 edges. 

The same d-periodic graph allows other placements. One may start with 
an arbitrary simplex PqPi...P (1 with the origin O in the interior. Then O is 
connected with the vertices Pi and provides d -I- 1 edge representatives. By 
using the periodicity lattice A generated as above by the edge vectors of the 
initial simplex, one extends the finite system of vertices and edges to a full 
A-periodic system. The corresponding frameworks will be called diamantine 
frameworks, frameworks of diamond type or simply diamond frameworks. 

The three dimensional standard framework corresponds with the bonded 
structure of carbon atoms in (ideal) diamond crystals. The two dimensional 
version may be imagined as a graphene layer. 
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What distinguishes the standard diamond framework among its diamantine 
relatives is its full symmetry group, which is maximal. We refer to cniEii 
[27] for the corresponding theory of canonical placements of periodic graphs. 
Related aspects can be found in nmnisi. Our present inquiry will assume 
that a diamantine framework in has been given and will be concerned with 
its deformations as a periodic bar-and-joint structure. Deformations will have 
new placements for the vertices (i.e. joints) of the structure, but in such a 
way that all edges preserve their length (i.e. behave like rigid bars) and the 
deformed framework remains periodic with respect to the abstract periodicity 
group marked at the outset [5j. It is important to retain the fact that the 
representation of this abstract periodicity group by a lattice of translations of 
rank d is allowed to vary in deformations. 


2 Deformations 

Let po, ■■■Pd be the edge vectors emanating from the origin. Their squared norm 
will be denoted by {pi,Pi) = Si. The lattice of periods is generated by pi — 
Pq, i = 1,d and with this ordering, the oriented volume of the fundamental 
parallelepiped (unit cell) is given by 


Vj(p) = det 


1 1 

Pio Pii 


1 

Pid 


det [ p,^ - Pig ... p,^ -p,g ] (1) 


The deformation space of a diamond framework defined by p, with V{p) 0 
can be parametrized by the connected component of p in 

d 

X{S^^]\{p-.V{p) = 0} 

i=0 

modulo the natural action of SO{d). As usual, denotes the (d — 1)- 

dimensional sphere of radius r. 

We investigate first the critical points of V{p) on the indicated product of 
spheres. At a critical point p, we must have 

dV 

— (p) = X,p^ 
api 

for some Lagrange multipliers A^, i = 0, ...,d. By inner product with pj, we 
obtain: 


and for i ^ j this gives: 



{P),P3) = K{Pi,Pj) 


^i{Pi,Pj) = det\po...pi...pd\ 


( 2 ) 

( 3 ) 
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where pi means the absence of that column. 
For i = j we find: 


WPi? = v{p) + {-iy^^det\pQ...pi...pd\ (4) 

By summation or by Euler’s formula for the homogeneous function V (p), we 
obtain: 


^ = dV{p) 


2=0 


(5) 


Equations (§ and show that for non-zero critical values, all Ai ^ 0 and then 
it follows that all inner products {pi,pj), i j must be equal. If we denote by 
a their common value, we see from the Gram determinant of the vectors pi that 
a must be a root of the degree (d + 1) equation 



( So 

a 

a 

. a 

det 

a 

Si 

a 

. a 


\ a 

a 

a . 

■ Sd 


( 6 ) 


For a = Si = \pi\^, the determinat on the left hand side of equals Si Y\j^i{sj — 
Si). If we assume distinct values Si and the natural ordering sq < Si < < Sd, 

we see that equation (§ must have a root in every interval (si, Si+i), since the 
determinant has opposite signs at the endpoints. Thus, in the general case of 
distinct Si, equation ([^ has d positive simple roots. The remaining root must 
be negative since the term of degree one in a is zero i.e. the sum of all products 
of d roots vanishes. By continuity, we still have a unique negative root a_ when 
some of the values Si become equal. 


With respect to positive roots of (§, we note that the inequality {pi,pj) < \pi\- 
\Pj \ — implies that only the smallest among them a+ may correspond 

to a cricical point. Moreover, from V{p) ^ 0, we infer that this root must be 
simple. Indeed, otherwise, with tso = Si < S 2 < ... < Sd, the root must be 
So = Si = (PoiPi) and this forces po = Pii resulting in V{p) = 0. In fact, for the 
simple root q:+ one may confirm the more precise location a+ S (so,sy^s}^^) 
by evaluating the determinant in ([^ at a = and finding it negative. 

With ri = and a = rgri, we have 


/ So 


det 


a Si a 

Y a a a 


a \ 
a 

Sd j 


( 1 
1 


= sosidet 


?’i 

’’o 


ri ro S2 a... 


V ri ro 


ri \ 
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/O 1 1 ... 1 \ 

= soSi(ro - ri)^dei 1 S2 a ... a ^ 

\ 1 a a ... Sd J 

= -soSi(ro - rif TT(si - a)(V —-—) < 0. 

Si — a 

i=2 i=2 * 

Observing that, up to SO{d), a Gram matrix detrmines two reflected vector 
configurations with corresponding volumes of opposite signs and equal absolute 
value, we may eatablish by induction on d > 2, that the negative root corre¬ 
sponds to the absolute maximum and minimum, while a simple smallest positive 
root corresponds to reflected critical configurations which are neither local min¬ 
ima nor local maxima. For d = 2, the deformation space is a surface and a 
‘seddle point’ configuration is illustrated in Figure Note that, in absolute 
value, a local variation of po increases the area of triangle poPiP2i while a local 
variatin of pi decreases the area. 



Figure 4: A configuration p = (j>o,Pi,P 2 ) corresponding to a saddle point. 
Note that the origin O is the orthocenter of triangle PoPiP 2 ■ 


Let us observe, for the latter part of the argument, when sq < si < S 2 < ... < Sd 
and {pi,Pj) = € (sqj i j, that in the hyperplane passing through 

Pi, ...,pd we obtain a (d— l)-dimensional configuration with d vectors qi, where 

Pi = p-Po+qi, (Po,di)=0, i = l,...,d (7) 
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with fipo in the indicated hyperplane. Since p, = a+/so, we find 

Oi-L 

{qr,qj) = a+{l -) < 0 

So 

q,2 

which is a maximum or minimum configuration for {qi,qi) = Si —> 0. 

The absence of local minima or maxima other than the absolute ones shows 
that 11^=0 \ {P • ^(p) “ exactly two connected components which 

are isomorphically exchanged by any orientation-reversing orthogonal transfor¬ 
mation. 

Theorem 1 Let Sq < Si <,... < denote the squared lengths of a complete 
set of edge representatives for a d-dimensional diamond framework 'D{p). The 
deformation space of this framework is a manifold of dimension — 1 and 

can be described as the quotient nto 5^^ \{p :V{p) = 0}/O{d). The squared 
volume of a fundamental parallelepiped (unit cell) has a maximum corresponding 
to the unique negative root a_ of &■ 

The indicated quotient space is a manifold because 0{d) acts without fixed 
points on V{p) 0 and local charts can be obtained as transversals to orbits. 
An alternative argument will follow from using Gram coordinates. 

We have already implicated in our considerations the image of the parameter 
space 11^=0*^^^ \ {P ■ ^{p) = ^}/0{d) obtained via Gram matrices G{p) = 
i(PijPj))o<i,j<d, with prescribed diagonal entries {pi,Pi) = Si, i = 0, ...,d. The 
image is the locus of all semipositive (d -I- 1) x (d -f 1) symmetric matrices of 
rank d with positive diagonal minors. 

While it is generally true that deformation spaces of periodic frameworks are 
connected components of semi-algebraic sets [6], the case of diamond frameworks 
allows very explicit treatment. We are going to describe another image of the 
deformation space, given by the Gram marices of the generators Vi = Pi — Po, 
i = 1 ,..., d of the corresponding periodicity lattices. 

We denote by w = the symmetric d x d matrix with entries 

ujij = {vi,Vj) = (p^-po,Pj -po) = {Pi,Pj) - {p^+PJ,Po) + So (8) 

From ujii = Si -\- sq — 2(pi,po) we retrieve all entries 

{Pi,Po) = ^{si + So - uJii), i = l,...,d (9) 

and then 

{Pi,Pj) = -\-^{si-\-Sj - OJU - ojjj), l<if=j<d (10) 
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Thus, G{p) and uj contain equivalent information and the ajfine transformation 
defined by ([^ and (10) take the algebraic hypersurface detG{p) = 0 to a degree 
(c? + 1) hypersurface in the vector space oi d x d symmetric matrices (with 
coordinates ujij,1 < i < j < d). Line and column operations give the equivalent 
expression as a ‘bordered determinant’ of w, namely 


/ 

det 

V 


So 

5(si — So — Wii) 
2 (Srf So ^dd) 


|(si — So — Wii) 
Wll 

OJdl 


2 (StZ So ^dd) \ 
Wdd / 


0 ( 11 ) 


The image of the deformation space consists of the intersection of this hyper¬ 
surface with the open cone of positive definite symmetric matrices. 


Remark: Strictly speaking, we should indicate by an index s = (so,...,Sd) 
the fact that the diagonal entries in the Gram matrix G{p) = Gs{p) are fixed 
by {pi,Pi) = Si- While all hypersurfaces detGs(p) = 0 in coordinates Oij = 
{piiPj)i 0 ^ * < j < d are isomorphic to the case Si = 1, i = 0,...,d under 
linear transformations of the form Oij i—aij/(siSj)^/^, what varies with s is 
the intersection with the locus of semipositive definite matrices with positive 
diagonal minors (or, in terms of coordinates Wij in (11), the intersection with 
the positive definite cone). 


Let us show here that for all s, Si > 0,i = 0,...,d, there are no singularities 
of detGa{p) = 0 in this intersection. With Oij = {pi,pj) as coordinates, the 
vanishing of the gradient would make all off-diagonal minors equal to zero. 
Since Gs{p) is not invertible, the diagonal minors must vanish as well and we 
obtain a contradiction. The argument shows in fact that the singularity locus 
of detGs ip) = 0 is made of symmetric matrices of rank strictly less than d. 


We shall investigate in more detail the case d = 2, which involves Cayley nodal 
cubic surfaces. 


3 Planar diamantine frameworks 

For d = 2, the deformation space of a diamantine framework is a surface. The 
topology of this surface can be determined with relative ease by observing that 
in the quotient 


f[Sl^JSO{2) 

i=0 

the action of «S'0(2) = can be used to fix p 2 . Thus, we may use as parameters 
the angles ( 00 , 4>i) G x of po and pi with p 2 and exclude from this torus, 
the degenerate locus corresponding to collinear vertices po,pi,p 2 . 






When So < Si < S 2 the torus is cut by two disjoint loops into two open cylin¬ 
ders (isomorphic by reflection in the p 2 axis), The generic case is illustrated in 
Figure 



Figure 5: Topology determination in the generic case sq < si < S 2 . With 
P 2 fixed, the possibilities for {po,Pi) are parametrized by the torus x S^. 
The picture indicates the projection on the secon factor parametrizing pi. The 
two dark loops depicted on the torus represent degenerate configurations with 
collinear vertices po,pi,p 2 ■ 

For So = Si < ls 2 we cut along intersecting loops which leave two (reflection 
isomorphic) topological open discs in the complement. The case s-|-0 = si =S 2 
is the simplest, since we exclude (po = 0, (l)i — 0 and 4>o = pi, leaving two 
topological open discs in the complemeent. This gives: 

Proposition 2 The deformation space of a planar diamantine framework is an 
open cylinder when sq < si < S 2 and an open disc when sq = si < S 2 . 

Remark: The topological change from disc to cylinder is reflected in the ap¬ 
parition of the saddle point for the area function. 

With these preparations, we may follow the scenario described in the previous 
section whereby the deformation space is mapped to the positive definite cone 
(of symmetric 2x2 matrices) by associating to a framework the Gram matrix 
u! = uj{p) of the two generators of the periodicity lattice Pi — Pq, i = 1,2. 

We illustrate the ‘standard’ case sq = si = S 2 = 1 which allows the most 
symmetric presentation, already depicted in Figure Equation (|TI|) takes the 


form 


fi^) — /(Wl 1)1^12,^22) — 



( 12 ) 
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Figure 6: The Cayley cubic surface and the positive semidefinite cone passing 
through its four nodes (rendered as small gaps between nonsingular connected 
components) . The 2-diamond deformation surface is the part inside the cone 
which wraps the tetrahedron with vertices at the nodes. The three lines of 
tangency along the three edges from the origin are excluded. 


The gradient vector 

( ^22(5^12 — 5W11 — ^ UJ22 -f 1) \ 

^UJiiOJ22 — 2uJi2 j (13) 

^11(5^12 — 5W22 — |wil + 1) / 

vanishes at the four points: (0,0,0), ((4,0,0), (0,0,4), (4,4,4). At this stage, 
one may recognize that our cubic (HI is projectively equivalent with the Cayley 
nodal cubic surface usually given in projective coordinates (xg : Xi : X 2 ■ xg) by 
the vanishing of the third symmetric polynomial 


a 3 {x) = X 1 X 2 X 3 + X 0 X 2 X 3 + X 0 X 1 X 3 + X 0 X 1 X 2 = 0. (14) 

Up to projective equivalence, , the Cayley cubic is completely characterized by 
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the fact that it has exactly four singularities which are all nodal points (which 
may be seen as the vertices of a tetrahedron whose edge-supporting lines belong 


to the cubic). Although equation (14) immediately reveals projective symmetry 
under the permutation group of the four coordinates Xi, we shall retain the 


description (12) in terms of w, since we have to intersect our cubic with the 


positive definite cone defined by 


tr{uj) = ojii + UJ 22 > 0 and det{uj) = (JJ 11 OJ 22 — ^^->12 > -0 


(15) 


The bounding cone det(ijj) = 0 intersects the cubic (12) along the three double 
lines given by the edges through (0, 0,0) of the four nodes tetrahedron. More¬ 
over, detipj) = 0 is also the tangent cone at the singularity (0,0,0) for (12). 


4 Auxetic trajectories 

The notion of auxetic behaviour is formulated in the context of elasticity theory 
and is an expression of negative Poisson’s ratios miiiin]. Simply phrased, 
auxetic behaviour means becoming laterally wider when streched and thinner 
when compressed. Our purely mathematical treatment of periodic framework 
deformations cannot implicate physical properties per se and we shall rely on 
the geometrical alternative proposed in [7] and illustrated in [5] for periodic 
pseudo-triangulations. 

The auxetic property, from the greek word for growth, refers to certain one- 
parameter deformations of a periodic framework and not to the framework itself, 
which may allow, in general, non-auxetic deformations as well. We use the 
following definition which involves the positive semidefinite cone f2((i) in the 
vector space oi d x d symmetric matrices i.e. the cone made of all symmetric 
matrices with non-negative eigenvalues. 

Definition 3 A smooth one-parameter deformation of a periodic framework in 
is called an auxetic path, or simply auxetic, when the corresponding curve of 
Gram matrices w(t) for a set of independent generators of the periodicity lattice 
has all its tangent vectors duj{T)/dT in the positive semidefinite cone of d x d 
symmetric matrices. 


Remark. The similarity with the notion of causal-line in special relativity is 
apparent. A causal-line in a Minkowski space is a smooth curve with all its 
tangents in the light cone. For d = 2, the boundary of the positive semidefinite 
cone Cl(2) is given by the vanishing of a single (Lorenzian) quadratic form and 
the two notions coincide. For d > 2, the geometry of the cone Cl(d) is more 
complicated [HE], but the analogy persists. 

We shall make use of the description obtained in Section for the deforma¬ 
tion surface of 2-diamnond and investigate the possibility of auxetic trajectories 
on this surface. 
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In order to identify the region of this surface made of points allowing some non¬ 
trivial auxetic trajectory through them, we must look for points where some 
non-zero tangent vector belongs to the positive semidefinite cone (as a free 
vector). The boundary of this region consists therefore of points with tangent 
planes parallel to a tangent plane of the cone surface along a generating line. 
This locus can be detected algebraically as follows. 

We use the natural Euclidean norm 11 11 1 j. defined on the space of 2 x 2 symmetric 
matrices by 


|w||tr = + 2a;j2 +<^ 22 )^^^ 


(16) 


The positive semidefinite cone is self-dual with respect to this norm m. hence 
verifying that a tangent plane corresponts to a tangent plane of the cone det{u}) = 
0 along a generating line amounts to verifying that the normal direction with 


respect to (161 satisfies the cone equation. The required normal is immediately 


obtained from the gradient formula (13) by halving the middle coordinate. Thus, 


the equation for the region’s boundary is the quartic 


, . ,111 ..111 
g{w) = a;iia; 22 (-a;i 2 - -wn - -0122 + l)( 2 '^i 2 - ^^^22 - ^wii + 1) 


— (-a;iia;22 — 1 ^ 12 )^ = 0 (17) 

The curve of degree twelve resulting from the intersection of the Cayley cubic 
and the above quartic can be first guessed and then proven to be made of the six 
lines supporting the edges of the tetrahedon of nodes, counted with multiplicity 
two. For instance, intersecting the quartic 0 with the plane UJ 12 = 0 gives the 
four lines resulting from the factorization 

UJ 22 ) = gWiiW 22 (wii -I- a ;22 — 4)(a;ii -I- UJ 22 — 2) (18) 

The first three factors correspond with the three edges of the tetrahedron of 
nodes on the face uji 2 . By symmetry and the fact that auxetic trajectories 
cannot exist in a neighborhood of the standard 2-diamond, we conclude that 
the region of the deformation surface consisting of points allowing a non-trivial 
auxetic deformation path is given by intersecting with the open half-space 


wii — UJ12 + 0J22 < 4 


(19) 


If we express (191 in terms of the unit vectors pi, j = 0,1, 2, we find the condition 


- 1 < (Po,Pi)-k (pi,P2) + (P2,Po) (20) 

Geometrically, this means pointedness: the three (edge) vectors pi aree all on 
one side of some line through the origin. Equivalently, one of the three sectors 


12 







Figure 7: Three unit vectors in a pointed configuration (A), boundary (tran¬ 
sition) case (B) and non-pointed configuration (C), corresponding to <,= and 
> in relation (20). 


determined around the origin by the three vectors has an angle larger than tt 
as illustrated on the left configuration of Figure 

We summarize our conclusion in the following statement. 

Proposition 4 A framework obtained from the standard 2-diamond framework 
allows some non-trivial (small) auxetic deformation path if and only if pointed 
at all vertices. 

Remarks. We have argued above for pointedness at the vertex of reference. 
Diamantine frameworks are vertex transitive under full crystallographic sym¬ 
metry, hence pointedness at all vertices follows. The reentrant honeycomb in 
Figure is pointed at all vertices. For the relevance of pointedness in periodic 
pseudo-triangulations we refer to [ilTlIH]. 
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